BaBar-PUB-11/014 

SLAC-PUB-14811 

hep-ex/1112.0702 



Study of 5 —7- Xuiiy decays in BB events tagged by a fully reconstructed S-meson 

decay and determination of \Vub\ 

J. P. Lees, V. Poireau, and V. Tisserand 
Laboratoire d'Annecy-le-Vieux de Physique des Particules (LAPP), 
Universite de Savoie, CNRS/IN2P3, F-74941 Annecy-Le-Vieux, France 

J. Garra Tieo and E. Granges 
Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 

M. Martinelli"^ D. A. Milanes", A. Palano''^ and M. Pappagallo''^ 
INFN Sezione di Barif^ ; Dipartimento di Fisica, Universita di Barf , 1-70126 Bari, Italy 

G. Eigen and B. Stugu 
University of Bergen, Institute of Physics, N-5007 Bergen, Norway 

D. N. Brown, L. T. Kerth, Yu. G. Kolomensky, G. Lynch, and K. Tackmann 
Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 

H. Koch and T. Schroeder 
Ruhr Universitat Bochum, Institut fiir Experimentalphysik 1, D-44780 Bochum, Germany 

D. J. Asgcirsson, C. Hearty, T. S. Mattison, and J. A. McKenna 
University of British Columbia, Vancouver, British Columbia, Canada V6T IZl 

A. Khan 

Brunei University, Uxbndge, Middlesex UBS 3PH, United Kingdom 

V. E. Bhnov, A. R. Buzykaev, V. P. Druzhinin, V. B. Golubev, E. A. Kravchenko, A. P. Onuchm, 
S. L Serednyakov, Yu. L Skovpen, E. P. Solodov, K. Yu. Todyshev, and A. N. Yushkov 
Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

M. Bondioh, D. Kirkby, A. J. Lankford, M. Mandelkern, and D. P. Stoker 
University of California at Irvine, Irvine, California 92697, USA 

H. Atmacan, J. W. Gary, F. Liu, O. Long, and G. M. Vitug 

University of California at Riverside, Riverside, California 92521, USA 

C. Campagnari, T. M. Hong, D. Kovalskyi, J. D. Richman, and C. A. West 
University of California at Santa Barbara, Santa Barbara, California 93106, USA 

A. M. Eisner, J. Kroseberg, W. S. Lockman, A. J. Martinez, T. Schalk, B. A. Schumm, and A. Seiden 

University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 

C. H. Cheng, D. A. Doll, B. Echenard, K. T. Flood, D. G. Hithn, P. Ongmongkolkul, F. C. Porter, and A. Y. Rakitin 

California Institute of Technology, Pasadena, California 91125, USA 

R. Andrcassen, M. S. Dubrovin, Z. Huard, B. T. Meadows, M. D. Sokoloff, and L. Sun 

University of Cincinnati, Cincinnati, Ohio 45221, USA 



P. C. Bloom, W. T. Ford, A. Gaz, M. Nagel, U. Nauenberg, J. G. Smith, and S. R. Wagner 
University of Colorado, Boulder, Colorado 80309, USA 



2 

R. Ayad* and W. H. Toki 

Colorado State University, Fort Collms, Colorado 80523, USA 

B. Spaan 

Techmsche Umversitdt Dortmund, Fakultdt Physik, D-44221 Dortmund, Germany 

M. J. Kobel, K. R. Schubert, and R. Schwierz 
Technische Umversitdt Dresden, Institut fiir Kern- und Teilchenphysik, D-01062 Dresden, Germany 

D. Bernard and M. Verderi 

Laboratoire Leprince-Ringuet, Ecole Polytechmque, CNRS/IN2P3, F-91128 Palaiseau, France 

P. J. Clark and S. Playfer 

University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 

D. Bettoni", C. Bozzi", R. Calabrese'^\ G. Cibinetto''\ E. Fioravanti"^ 1. Garzia"^ 
E. Luppi°^ M. Munerato"^ M. Negrini"^ A. Petrella°^ L. Piemontese", and V. Santoro 
INFN Sezione di Ferrara"^ ; Dipartimento di Fisica, Universitd di Ferrard' , 1-44100 Ferrara, Italy 

R. Baldini-Ferroli, A. Calcaterra, R. de Sangro, G. Finocchiaro, 
M. Nieolaei, P. Patteri, I. M. Peruzzi,'!' M. Piccolo, M. Rama, and A. Zallo 
INFN Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 

R. Contri"^ E. Guido''^ M. Lo Vetere''^ M. R. Monge''^ S. Passaggio", C. Patrignani''^ and E. Robutti'' 
INFN Sezione di Genova"^ ; Dipartimento di Fisica, Universitd di Genoval' , I-I6I46 Genova, Italy 

B. Bhuyan and V. Prasad 
Indian Institute of Technology Guwahati, Guwahati, Assam, 781 039, India 

C. L. Lee and M. Morii 

Harvard University, Cambridge, Massachusetts 02138, USA 

A. J. Edwards 
Harvey Mudd College, Claremont, California 91711 

A. Adametz, J. Marks, and U. Uwer 
Universitat Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 

F. U. Bernlochner, M. Ebert, H. M. Lacker, and T. Lueck 
Humboldt-Universitdt zu Berlin, Institut fiir Physik, Newtonstr, 15, D-12489 Berlin, Germany 

P. D. Dauncey and M. Tibbetts 
Imperial College London, London, SW7 2AZ, United Kingdom 

P. K. Behera and U. Mallik 
University of Iowa, Iowa City, Iowa 52242, USA 

C. Chen, J. Cochran, W. T. Meyer, S. PreU, E. L Rosenberg, and A. E. Rubin 

Iowa State University, Ames, Iowa 5001 1-3160, USA 

A. V. Gritsan and Z. J. Guo 
Johns Hopkins University, Baltimore, Maryland 21218, USA 

N. Arnaud, M. Davier, G. Grosdidier, F. Le Diberder, A. M. Lutz, 
B. Malaescu, P. Rondeau, M. H. Schune, A. Stocchi, and G. Wormser 
Laboratoire de lAccelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, 
Centre Scientifique d'Orsay, B. P. 34, F-91898 Orsay Cedex, France 



3 



D. J. Lange and D. M. Wright 
Lawrence Livermore National Laboratory, Livermore, California 94550, USA 

I. Bingham, C. A. Chavez, J. P. Coleman, J. R. Fry, E. Gabathuler, D. E. Hutchcroft, D. J. Payne, and C. Touramanis 

University of Liverpool, Liverpool L69 7ZE, United Kingdom 

A. J. Bevan, F. Di Lodovico, R. Sacco, and M. Sigamani 
Queen Mary, University of London, London, El 4NS, United Kingdom 

G. Cowan 

University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 

D. N. Brown and C. L. Davis 
University of Louisville, Louisville, Kentucky 40292, USA 

A. G. Dcnig, M. Fritsch, W. Gradl, A. Hafner, and E. Prcncipc 
Johannes Cutenberg-Umversitdt Mamz, Institut fiir Kernphysik, D-55099 Mainz, Cermany 

K. E. Alwyn, D. Bailey, R. J. Barlow, G. Jackson, and G. D. Lafferty 
University of Manchester, Manchester MIS 9PL, United Kingdom 

R. Ccnci, B. Hamilton, A. Jawahcry, D. A. Roberts, and G. Simi 
University of Maryland, College Park, Maryland 20742, USA 

C. Dallapiccola 

University of Massachusetts, Amherst, Massachusetts 01003, USA 

R. Cowan, D. Dujmic, and G. ScioUa 

Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 

D. Lindcmann, P. M. Patel, S. H. Robertson, and M. Schram 

McCill University, Montreal, Quebec, Canada H3A 2T8 

P. Biassoni''^ A. Lazzaro'^\ V. Lombardo'', N. Ncri'^^ F. Palombo''^ and S. Stracka"* 
INFN Sezione di Milano"' ; Dipartimento di Fisica, Universitd di Milano^ , 1-20133 Milano, Italy 

L. Cremaldi, R. Godang,^ R. Krocger, P. Sonnek, and D. J. Summers 
University of Mississippi, University, Mississippi 38677, USA 

X. Nguyen and P. Taras 
Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada II3C 3J7 

G. De Nardo''^ D. Monorchio'^^ G. Onorato"^ and C. Sciacca'^'' 
INFN Sezione di Napoli^ ; Dipartimento di Scienze Fisiche, 
Universitd di Napoli Federico if, 1-80126 Napoli, Italy 

G. Raven and H. L. Snoek 
NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

C. P. Jessop, K. J. Knoepfel, J. M. LoSecco, and W. F. Wang 
University of Notre Dame, Notre Dame, Indiana 46556, USA 

K. Honschcid and R. Kass 
Ohio State University, Columbus, Ohio 43210, USA 

J. Brau, R. Frey, N. B. Sinev, D. Strom, and E. Torrence 
University of Oregon, Eugene, Oregon 97403, USA 



4 



E. Feltresi°^ N. Gagliardi°^ M. Margoni°^ M. Morandin", 
M. Posocco'", M. Rotondo", F. Simonctto^^ and R. Stroili"'' 
INFN Sezione di Padova^ ; Dipartimento di Fisica, Universitd di Padovci , 1-35131 Padova, Italy 

E. Ben-Haim, M. Boniben, G. R. Bonneaud, H. Briand, G. Caldcrini, 
J. Chauveau, O. Hamon, Ph. Leruste, G. Marchiori, J. Ocariz, and S. Sitt 
Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, Universite Pierre et Marie Curie-Paris6, 
Universite Denis Diderot-Paris?, F-75252 Pans, France 

M. Biasini'^^ E. Manoni'^^ S. Pacetti'^^ and A. Rossi"'' 
INFN Sezione di Perugia'^ ; Dipartimento di Fisica, Universita di Perugia'' , 1-06100 Perugia, Italy 

C. Angelini°^ G. Batigna^i'^^ S. Bettarini'^^ M. Cal•pinelli°^^ G. Casarosa°^ A. Cervelli"^ F. Forti'^^ 
M. A. Giorgi"^ A. Lusiani'^^ B. ObcrhoP^ E. Paolo^i°^ A. Perez", G. Rizzo°^ and J. J. Walsh" 
INFN Sezione di Pisa"^ ; Dipartimento di Fisica, Universitd di Pisa'' ; Scuola Normale Superiore di Pisa'^, 1-56127 Pisa, Italy 

D. Lopes Pegna, C. Lu, J. Olsen, A. J. S. Smith, and A. V. Tehiov 

Princeton University, Princeton, New Jersey 08544; USA 

F. AnuUi", G. Gavoto'', R. Faccini"^ F. Ferrarotto", F. Ferroni"^ 

M. Gaspero''^ L. Li Gioi°, M. A. Mazzoni", and G. Piredda" 
INFN Sezione di Roma"; Dipartimento di Fisica, 
Universita di Roma La Sapienza'' , 1-00185 Roma, Italy 

G. Biinger, O. Griinberg, T. Hartmann, T. Leddig, H. Schroder, and R. Waldi 
Universitdt Rostock, D- 18051 Rostock, Germany 

T. Adye, E. O. Olaiya, and F. F. Wilson 

Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 

S. Emery, G. Hamel de Monchenault, G. Vasseur, and Gh. Yeche 
CEA, Irfu, SPP, Centre de Saclay, F-91191 Gif-sur-Yvette, France 

D. Aston, D. J. Bard, R. Bartoldus, G. Gartaro, M. R. Gonvery, J. Dorfan, G. P. Dubois-Felsmann, W. Dunwoodie, 
R. G. Field, M. Franco Seviha, B. G. Fulsom, A. M. Gabareen, M. T. Graham, P. Grenier, G. Hast, W. R. Innes, 
M. H. Kelscy, H. Kim, P. Kim, M. L. Kocian, D. W. G. S. Leith, P. Lewis, S. Li, B. Lindquist, S. Luitz, V. Luth, 
H. L. Lynch, D. B. MacFarlane, D. R. MuUer, H. Neal, S. Nelson, L Ofte, M. Perl, T. PuUiam, B. N. Ratcliff, 
A. Roodman, A. A. Salnikov, R. H. Schindler, A. Snyder, D. Su, M. K. Suhivan, J. Va'vra, A. P. Wagner, 
M. Weaver, W. J. Wisniewski, M. Wittgen, D. H. Wright, H. W. Wulsin, A. K. Yarritu, G. G. Young, and V. Ziegler 
SLAC National Accelerator Laboratory, Stanford, California 94309 USA 

W. Park, M. V. Purohit, R. M. White, and J. R. Wilson 
University of South Carolina, Columbia, South Carolina 29208, USA 

A. Randle-Gonde and S. J. Sekula 

Southern Methodist University, Dallas, Texas 75275, USA 

M. Belhs, J. F. Benitez, P. R. Burchat, and T. S. Miyashita 
Stanford University, Stanford, California 94305-4060, USA 

M. S. Alam and J. A. Ernst 

State University of New York, Albany, New York 12222, USA 

R. Gorodeisky, N. Guttman, D. R. Peimer, and A. Soffer 
Tel Aviv University, School of Physics and Astronomy, Tel Aviv, 69978, Israel 



P. Lund and S. M. Spanier 



5 



University of Tennessee, Knoxville, Tennessee 37996, USA 

R. Eckmann, J. L. Ritchie, A. M. Ruland, C. J. Schilling, R. F. Schwittcrs, and B. C. Wray 
University of Texas at Austin, Austin, Texas 78712, USA 

J. M. Izen and X. C. Lou 

University of Texas at Dallas, Richardson, Texas 75083, USA 

F. Bianchi"'' and D. Gamba"'' 
INFN Sezione di Torino"" ; Dipartimento di Fisica Sperimentale, Universitd di Torino'' , 1-10125 Torino, Italy 

L. Lanceri'^'' and L. Vitale"^ 
INFN Sezione di Trieste" ; Dipartimento di Fisica, Universitd di Trieste'' , 1-34127 Trieste, Italy 

V. Azzolini, F. Martinez- Vidal, and A. Oyanguren 
IFIC, Universitat de Valencia-CSIC, E-46071 Valencia, Spain 

H. Ahmed, J. Albert, Sw. Banerjee, H. H. F. Choi, G. J. King, R. Kowalewski, 
M. J. Lewczuk, C. Lindsay, L M. Nugent, J. M. Roney, R. J. Sobie, and N. Tasneem 

University of Victoria, Victoria, British Columbia, Canada V8W 3P6 

T. J. Gershon, P. F. Harrison, T. E. Latham, and E. M. T. Puccio 

Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 

H. R. Band, S. Dasu, Y. Pan, R. Prepost, and S. L. Wu 
University of Wisconsin, Madison, Wisconsin 53706, USA 

We report measurements of partial branching fractions for inclusive charmless semileptonic B 
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e'^e~ storage rings. We select events in which the decay of one of the B mesons is fully reconstructed 
and an electron or a muon signals the semileptonic decay of the other B meson. We measure partial 
branching fractions AB in several restricted regions of phase space and determine the CKM element 
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I. INTRODUCTION 

A principal physics goal of the BaBAR experiment is 
to establish CP violation in B meson decays and to test 
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whether the observed effects are consistent with the Stan- 
dard Model (SM) expectations. In the SM, CP-violating 
effects result from an irreducible phase in the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix [1, 2]. 
Precise determinations of the magnitude of the matrix 
element \Vub\ will permit more stringent tests of the SM 
mechanism for CP violation. This is best illustrated in 
terms of the unitarity triangle [3], the graphical repre- 
sentation of one of the unitarity conditions of the CKM 
matrix, for which the side opposite to the angle /3 is pro- 
portional to the ratio The best way to deter- 
mine \Vub\ is to measure the decay rate for B Xuiv 
(here X refers to a hadronic final state and the index c 
or u indicates whether this state carries charm or not), 
which is proportional to |Ki5p. 
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There are two approaches to these measurements, 
based on either inclusive or exclusive measurements of 
semileptonic decays. The experimental uncertainties on 
the methods are largely independent, and the extraction 
of \Vub\ from the measured branching fractions relies on 
different sets of calculations of the hadronic contributions 
to the matrix element. For quite some time, the results 
of measurements of | Vub \ from inclusive and exclusive de- 
cays have been only marginally consistent [4, 5]. Global 
fits [6, 7] testing the compatibility of the measured an- 
gles and sides with the unitarity triangle of the CKM 
matrix reveal small differences that might indicate po- 
tential deviations from SM expectations. Therefore, it is 
important to perform redundant and improved measure- 
ments, employing different experimental techniques and 
a variety of theoretical calculations, to better assess the 
accuracy of the theoretical and experimental uncertain- 
ties. 

Although inclusive branching fractions exceed those of 
individual exclusive decays by an order of magnitude, the 
most challenging task for inclusive measurements is the 
discrimination between the rare charmless signal and the 
much more abundant decays involving charmed mesons. 
To improve the signal-to-background ratio, the events are 
restricted to selected regions of phase space. Unfortu- 
nately these restrictions lead to difficulties in calculat- 
ing partial branching fractions. They impact the conver- 
gence of Heavy Quark Expansions (HQE) [8, 9], enhance 
pcrturbativc and nonpcrturbative QCD corrections, and 
thus lead to significantly larger theoretical uncertainties 
in the determination of | Vub \ ■ 

We report herein measurements of partial branching 
fractions (AS) for inclusive charmless semileptonic B 
meson decays, B XuiE' [10]. This analysis extends 
the event selection and methods employed previously by 
BaBAR to a larger dataset [11]. We tag r{4S) BB 
events with a fully reconstructed hadronic decay of one 
of the B mesons (B^cco)- This technique results in a low 
event selection efficiency, but it uniquely determines the 
momentum and charge of both B mesons in the event, re- 
ducing backgrounds significantly. For charged B mesons 
it also determines their flavor. The semileptonic decay 
of the second B meson (i?rocoii) is identified by the pres- 
ence of an electron or a muon and its kinematics are 
constrained such that the undetectable neutrino can be 
identified from the missing momentum and energy of the 
rest of the event. However, undetected and poorly re- 
constructed charged particles or photons lead to large 
backgrounds from the dominant B — )■ Xciv decays, and 
they distort the kinematics, e.g., the hadronic mass Mx 
and the leptonic mass squared q^. 

For the Broco sample, the two dominant background 
sources are non-BB events from continuum processes, 
e+e^ — > qq{j) with q = it, d, s, or c, and com- 
binatorial BB background. The sum of these two 
backgrounds is estimated from the distribution of the 
beam energy-substituted mass m^s, which takes the 
following form in the laboratory frame: rriEs — 



+ pB ■ Pbeams)^/£'^eams ^ Pb' Hcrc pB refers to 
the momentum of the -Broco candidate derived from the 
measured momenta of its decay products, Pbcams = 
(£'bcams,Pboains) to the four-momcutum of the colliding 
beam particles, and y/s to the total energy in the T(4S') 
frame. For correctly reconstructed iJroco decays, the dis- 
tribution peaks at the B meson mass, and the width of 
the peak is determined by the energy spread of the col- 
liding beams. The size of the underlying background is 
determined from a fit to the jties distribution. 

We minimize experimental systematic uncertainties, 
by measuring the yield for selected charmless semilep- 
tonic decays relative to the total yield of semileptonic 
decays B — >■ X£p, after subtracting combinatorial back- 
grounds of the -Broco selection from both samples. 

In order to reduce the overall uncertainties, measure- 
ment of the signal B — > XuiP decays is restricted to 
regions of phase space where the background from the 
dominant B — > Xcii^ decays is suppressed and theo- 
retical uncertainties can be reliably assessed. Specifi- 
cally, signal events tend to have higher charged lepton 
momenta in the _B-mcson rest frame {p}), lower Mx, 
higher q^, and smaller values of the light-cone momen- 
tum P+ = Ex — \px 1 7 where Ex and px are energy and 
momentum of the hadronic system X in the B meson 
rest frame. 

The observation of charged leptons with momenta ex- 
ceeding the kinematic limit for B — > Xcii' presented first 
evidence for charmless semileptonic decays. This was 
followed by a series of measurements close to this kine- 
matic limit [12-16]. Although the signal-to-background 
ratio for this small region of phase space is favorable, the 
theoretical uncertainties are large and difficult to quan- 
tify. Since then, efforts have been made to select larger 
phase space regions, thereby reducing the theoretical un- 
certainties. The Belle Collaboration has recently pub- 
lished an analysis that covers about 88% of the signal 
phase space [17], similarly to one of the studies detailed 
in this article. 

We extract \Vub\ from the partial branching fractions 
relying on four different QCD calculations of the par- 
tial decay rate in several phase space regions: BLNP 
by Bosch, Lange, Neubert, and Paz [18-20]; DGE, the 
dressed gluon exponentiation by Andersen and Gardi [21, 
22]; ADFR by Aglietti, Di Lodovico, Ferrara, and Riccia- 
rdi [23, 24]; and GGOU by Gambino, Giordano, Ossola 
and Uraltsev [25]. These calculations differ significantly 
in their treatment of perturbative corrections and the 
parameterization of nonpcrturbative effects that become 
important for the different restrictions in phase space. 

This measurement of \Vub\ is based on combined sam- 
ples of charged and neutral B mesons. In addition, we 
present measurements of the partial decays rates for B^ 
and B~ decays separately. The observed rates are found 
to be equal within uncertainties. We use this observation 
to set a limit on weak annihilation (WA), the process 
bu ^ £~ Di, which is not included in the QCD calculation 
of the B X£D decay rates. Since final state hadrons 
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originate from soft gluon emission, WA is expected to 
contribute to the decay rate at large values of [26-29] . 

The outline of this paper is as follows: a brief overview 
of the BaBAR detector, particle reconstruction and the 
data and Monte Carlo (MC) samples is given in Sec- 
tion II, followed in Section III by a description of the 
event reconstruction and selection of the two event sam- 
ples, the charmless semileptonic signal sample and the 
inclusive semileptonic sample that serves as normaliza- 
tion. The measurement of the partial branching fractions 
and their systematic uncertainties are presented in Sec- 
tions IV and V. The extraction of | Vub \ based on four sets 
of QCD calculations for seven selected regions of phase 
space is presented in Section VI, followed by the conclu- 
sions in Section VII. 



II. DATA SAMPLE, DETECTOR, AND 
SIMULATION 

A. Data sample 

The data used in this analysis were recorded with the 
BaBAR detector at the PEP-II asymmetric energy e+e" 
collider operating at the T{AS) resonance. The total 
data sample, corresponding to an integrated luminosity 
of 426 fb"^ and containing 467 million T{AS) BB 
events, was analyzed. 



B. The BABAR detector 

The BaBAR detector and the general event reconstruc- 
tion are described in detail elsewhere [30, 31]. For 
this analysis, the most important detector features are 
the charged-particle tracking, photon reconstruction, and 
particle identification. The momenta and angles of 
charged particles are measured in a tracking system con- 
sisting of a five-layer silicon vertex tracker (SVT) and a 
40-layer, small-cell drift chamber (DCH). Charged parti- 
cles of different masses are distinguished by their ioniza- 
tion energy loss in the tracking devices and by the DIRC, 
a ring-imaging detector of internally reflected Cherenkov 
radiation. A finely segmented electromagnetic calorime- 
ter (EMC) consisting of 6580 CsI(Tl) crystals measures 
the energy and position of showers generated by electrons 
and photons. The EMC is surrounded by a thin super- 
conducting solenoid providing a 1.5 T magnetic field and 
by a steel flux return with a hexagonal barrel section 
and two endcaps. The segmented flux return (IFR) is in- 
strumented with multiple layers of resistive plate cham- 
bers (RPC) and limited streamer tubes (LST) to identify 
muons and to a lesser degree K^. 



C. Single particle reconstruction 

In order to reject misidentified and background tracks 
that do not originate from the interaction point, we re- 
quire the radial and longitudinal impact parameters to 
be ro < 1.5 cm and |zo| < 10 cm. For secondary tracks 
from Kg — >■ 7r"'"7r~ decays, no restrictions on the impact 
parameter are imposed. The efficiency for the reconstruc- 
tion of charged particles inside the fiducial volume for 
SVT, DCH, and EMC, defined by the polar angle in the 
laboratory frame, 0.410 < 6'iab < 2.54 rad, exceeds 96% 
and is well reproduced by Monte Carlo (MC) simulation. 

Electromagnetic showers are detected in the EMC as 
clusters of energy depositions. Photons are required not 
to be matched to a charged track extrapolated to the po- 
sition of the shower maximum in the EMC. To suppress 
photons from beam-related background, we only retain 
photons with energies larger than 50 MeV. Clusters cre- 
ated by neutral hadrons {K^ or neutrons) interacting in 
the EMC are distinguished from photons by their shower 
shape. 

Electrons are primarily separated from charged 
hadrons on the basis of the ratio of the energy deposited 
in the EMC to the track momentum. This quantity 
should be close to 1 for electrons since they deposit 
all their energy in the calorimeter. Most other charged 
tracks are minimum ionizing, unless they shower in the 
EMC crystals. 

Muons are identified by a neural network that com- 
bines information from the IFR with the measured track 
momentum and the energy deposition in the EMC. 

The average electron efficiency for laboratory momenta 
above 0.5 GeV is 93%, largely independent of momen- 
tum. The average hadron misidentification rate is less 
than 0.2%. Within the polar-angle acceptance, the av- 
erage muon efficiency rises with laboratory momentum 
and reaches a plateau of about 70% above 1.4 GeV. The 
muon efficiency varies between 50% and 80% as a func- 
tion of the polar angle. The average hadron misidenti- 
fication rate is about 1.5%, varying by about 0.5% as a 
function of momentum and polar angle. 

Charged kaons are selected on the basis of information 
from the DIRC, DCH, and SVT. The efficiency is higher 
than 80% over most of the momentum range and varies 
with the polar angle. The probability of a pion to be 
misidentified as a kaon is close to 2%, varying by about 
1% as a function of momentum and polar angle. 

Neutral pious are reconstructed from pairs of photon 
candidates that are detected in the EMC and arc as- 
sumed to originate from the primary vertex. Photon 
pairs having an invariant mass within 17.5 MeV (corre- 
sponding to 2.5 a) of the nominal tt*^ mass are considered 
tt" candidates. The overall detection efficiency, including 
solid angle restrictions, varies between 55% and 65% for 
7r° energies in the range of 0.2 to 2.5 GeV. 

— > TT+TT^ decays are reconstructed as pairs of tracks 
of opposite charge with a common vertex displaced from 
the interaction point. The invariant mass of the pair is 
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required to be in the range 490 < 771^+,^- < 505 MeV. 



D. Monte Carlo simulation 

We use MC techniques to simulate the response of the 
BaBAR detector [32] and the particle production and de- 
cays [33] , to optimize selection criteria and to determine 
signal efficiencies and background distributions. The 
agreement of the simulated distributions with those in 
data has been verified with control samples, as shown 
in Section IV D; the impact of the inaccuracies of the 
simulation is estimated in Section V. 

The size of the simulated sample of generic BE events 
exceeds the BB data sample by about a factor of three. 
This sample includes the common B Xclv decays. 
MC samples for inclusive and exclusive B — X^t-v decays 
exceed the size of the data samples by factors of 15 or 
more. 

Charmless semileptonic B — > X^iv decays are sim- 
ulated as a combination of resonant three-body decays 
with Xu = TT, 77, 77', p, uj, and decays to nonresonant 
hadronic final states X„. The branching ratios assumed 
for the various resonant decays are detailed in Table I. 
Exclusive charmless semileptonic decays are simulated 
using a number of different parameterizations: for B — > 
Tviv decays we use a single-pole ansatz [34] for the de- 
pendence of the form factor with a single parameter mea- 
sured by BABAR [35]; for decays to pseudoscalar mesons 
•q and 77' and vector mesons p and w we use form fac- 
tor parameterizations based on light-cone sum calcula- 
tions [36, 37]. 

The simulation of the inclusive charmless semileptonic 
B decays to hadronic states with masses larger than 
2m^ is based on a prescription by De Fazio and Neu- 
bert (DFN) [38] for the triple-differential decay rate, 
St I dq^ dEg dsH {Ei refers to the energy of the charged 
lepton and sh = M^) with QCD corrections up to 
0{as). The motion of the b quark inside the B me- 
son is incorporated in the DFN formalism by convolving 
the parton-level triple-differential decay rate with a non- 
perturbative shape function (SF). This SF describes the 
distribution of the momentum fc_|_ of the b quark inside 
the B meson. The two free parameters of the SF are A^^ 
and Ai^^. The first relates the B meson mass to the 
6 quark mass, mf^ = ruB — vl®^, and Ai^^ is the average 
momentum squared of the b quark. The SF parameter- 
ization is of the form F{k^) = N{1 — xj^e^^^"'^, where 
X = k+/A^^ < 1 and a = -3(iSF)2/AiSF _ ^ r^^^ g^.^^ 
three moments Ai of the SF must satisfy the following 
relations: Aq = 1, Ai = and A2 = -Ai^^/3. 

The nonresonant hadronic state X„ is simulated with 
a continuous invariant mass spectrum according to the 
DFN prescription. The fragmentation of the X„ system 
into final state hadrons is performed by jetset [39] . The 
resonant and nonresonant components are combined such 
that the sum of their branching fractions is equal to the 
measured branching fraction for inclusive B — >■ X^tv de- 



cays [40] , and the spectra agree with the DFN prediction. 
In order to obtain predictions for different values of A^^ 
and Ai^^, the generated events are reweighted. 

TABLE I: Branching fractions and their uncertainties [40] for 
exclusive B — >■ Xulv decays 



mode B{B° ~> Xuiu) BjE" Xulu) 
B^-kO? (136 ± 7) • 10"'' (77 ± 12) • 10-*^ 
B ^ rjii? (64 ± 20) • 10"*^ 

B ^ plv (247 ± 33) • 10"^ (128 ± 18) • 10"' 
B^ujIu (115 ± 17) • 10"' 

B ^ 7]'lv (17 ±22) ■ lO"'' 



We estimate the shape of background distributions 
by using simulations of the process e^e" — !■ T{4S) 
BB with the B mesons decaying according to measured 
branching fractions [40]. 

For the simulation of the dominant background from 
B Xctv decays, we have chosen a variety of differ- 
ent form factor parameterizations. For B — )• D(.v and 
B D*£v decays we use parameterizations [41] based 
on heavy quark effective theory (HQET) [42-45]. In 
the limit of negligible charged lepton masses, decays to 
pseudoscalar mesons are described by a single form fac- 
tor for which the dependence is expressed in terms 
of a slope parameter We use the world average 

p\, = 1.19 ± 0.06 [46], updated with recent precise mea- 
surements by the BABAR Collaboration [47, 48]. De- 
cays to vector mesons are described by three form fac- 
tors, of which the axial vector form factor dominates. 
In the limit of heavy quark symmetry, their q^ depen- 
dence can be described by three parameters for which 
we use the most precise BABAR measurements [47, 49]: 
pI,^ = 1.20 ± 0.04 [47, 49], Ri = 1.429 ± 0.074, and 
i?2 = 0.827 ± 0.044 [49]. For the simulation of semilep- 
tonic decays to the four L=l charm states, commonly re- 
ferred to as D** resonances, we use calculations of form 
factors by Leibovich, Ligeti, Stewart, and Wise [50]. We 
have adopted the prescription by Goity and Roberts [51] 
for nonresonant B — >■ D^*^X£D decays. 



III. EVENT RECONSTRUCTION AND SIGNAL 
EXTRACTION 

A. Reconstruction of hadronic B decays tagging 
BB events 

T{4:S) — )■ BB events are tagged by the hadronic de- 
cays of one of the B mesons based on a semi-exclusive 
algorithm that was employed in an earlier analysis [11]. 
We look for decays of the type Breco D'^*^Y^, where 
D^*^ is a charmed meson (D°, £>+, D*^, or D*^) and 
y is a charged state decaying to up to five charged 
hadrons, pious or kaons, plus up to two neutral mesons 
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or tt'^). The following decay modes of D mesons 
are reconstructed: K~tt^ , if~7r+7r'', if ~7r+7r~7r+, 

KItt+-k- and D+ K-tt+tt+, A'-tt+tt+tt", A'!?7r+, 
K°TT+Tr+Tr-, K°TT+TT° with K° tt+tt". D* mesons 
are identified by their decays, ^ Z)'^7r+, D+tt", 

and I?*'' ^ D^tt'', £'°7. Pions and photons from D* de- 
cays are of low energy and therefore the mass difference 
AM = m^Dir) — m{D) serves as an excellent discrimina- 
tor for these decays. 

Of the 1113 i?reco decay chains that we consider, we 
only retain the 342 ones with a signal purity V = 
S/{S+B) > 20%, where S and B, derived from MC sam- 
ples, denote the signal and background yields. The kine- 
matic consistency of the i?ioco candidates with B meson 
decays is checked using toes and the energy difference, 
AE = {Pb ■ Pboams - s/2)/y^. Wc rcstrict the B,.cco 
mass to Toes > 5.22 GeV and require AE = GeV 
within approximately three standard deviations, where 
the AE resolution depends on the decay chain. If an 
event contains more than one Bj-eco candidate, the decay 
chain with the highest probability is chosen. For this 
purpose we define 

2 2 /^^^D<" -^'^DoV fAEY 
Xtotal = Xvortox + + • (1) 

Here the first term is taken from a vertex fit for tracks 
from i?reco decays, the second relates reconstructed and 
nominal masses [40], M^(.) and Af £,(.), of the charm 

mesons {D^ , D^,D*° or D*"^), with the resolution cr^(.) , 
and the third term checks the energy balance AE com- 
pared to its resolution ctab- The number of degrees of 
freedom is therefore defined as N"^"^ = N$°{^^ + 2. The 
resulting overall tagging efficiency is 0.3% for S^S" and 
0.5% for B+B- events. 



B. Selection of inclusive B — > X£u decays 

In order to minimize systematic uncertainties, we mea- 
sure the yield of selected charmless semileptonic decays in 
a specific kinematic region normalized to the total yield 
of semileptonic B — >■ X£D decays. Both semileptonic 
decays, the charmless and the normalization modes, are 
identified by at least one charged lepton in events that are 
tagged by a Srcco decay. Both samples are background- 
subtracted and corrected for efficiency. Using this nor- 
malization, the systematic uncertainties on the -Brcco re- 
construction and the charged lepton detection cancel in 
the ratio or are eliminated to a large degree. 

The selection criteria for the charmless and the total 
semileptonic samples are chosen to minimize the statisti- 
cal uncertainty of the measurement as estimated from a 
sample of fully simulated MC events that includes both 
signal and background processes. 

A restriction on the momentum of the electron or 
muon is applied to suppress backgrounds from secondary 



charm or r decays, photon conversions and misidenti- 
fied hadrons. This is applied to p^, the lepton momen- 
tum in the rest frame of the recoiling B meson, which is 
accessible since the momenta of the T(4iS') and the recon- 
structed B are known. This transformation is important 
because theoretical calculations refer to variables that 
arc Lorcntz-invariant or measured in the rest frame of 
the decaying B meson. We require to be greater than 
1 GeV, for which about 90% of the signal is retained. 

For electrons and muons the angular acceptance is de- 
fined as 0.450 < 6 < 2.473 rad, where 9 refers to the po- 
lar angle relative to the electron beam in the laboratory 
frame. This requirement excludes regions where charged 
particle tracking and identification are not efficient. We 
suppress muons from J/ip decays by rejecting the event 
if a muon candidate paired with any other charged track 
of opposite charge (and not part of -Brcco) results in an 
invariant mass of the pair that is consistent with the J ftp 
mass. A similar requirement is not imposed on electron 
candidates, because of the poor resolution of the corre- 
sponding J/ip peak. 

We also reject events if the electron candidate paired 
with any other charged track of opposite charge is con- 
sistent with a 7 — > e+e" conversion. 

A variety of processes contributes to the inclusive 
semileptonic event samples, i.e. candidates selected by 
a Srcco decay and the presence of a high momentum lep- 
ton. In addition to true semileptonic decays tagged by a 
correctly reconstructed -Bj-oco, we consider the following 
classes of backgrounds: 

• Combinatorial background: the Brcco is not cor- 
rectly reconstructed. This background originates 
from BB or continuum e^e~ — ?> qqij) events. In 
order to subtract this background, the yield of true 
i?reco decays is determined from an unbinned max- 
imum likelihood fit to the toes distribution (Sec- 
tion HID). 

• Cascade background: the lepton does not origi- 
nate from a semileptonic B decay, but from sec- 
ondary decays, for instance from D mesons, includ- 
ing Dg — > TV, or residual J/rp background. 

• T background: electrons or muons originate from 
prompt T leptons, primarily from B — > XtD decays. 

• Fake leptons: hadrons are misidcntified as leptons, 
primarily muons. 

The last three sources of background are combined and 
in the following are referred to as "other" background. 

C. Selection of inclusive B — )■ Xuii' decays 

A large fraction of i? — > Xcii> decays is expected to 
have a second lepton from cascade decays of the charm 
particles. In contrast, in B Xy^lv decays secondary 
leptons are very rare. Therefore, we enhance signal 
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events by selecting events with only one charged lepton 
having > 1 GeV. 

In semileptonic B meson decays, the charge of the pri- 
mary lepton is equal to the sign of the charge of the b 
quark. Thus for B^B~ events in which the -Brcco and 
the lepton originate from different B decays in the event, 
we impose the requirement QbQe < 0, where Qb is the 
charge of the b quark of the i?reco and Qi is the charge 
of the lepton. For B^B° events this condition does not 
strictly hold because of flavor mixing. Thus, to avoid a 
loss in efficiency, this requirement is not imposed. The 
hadronic state Xu in charmless semileptonic decays is re- 
constructed from all particles that are not associated with 
the -Breco candidate or the charged lepton. The measured 
four-momentum Px is defined as 

1 = 1 2 = 1 

where the summation extends over the four-vectors of 
the charged particles and photon candidates. From this 
four-vector, other kinematic variables, = = 

- Pxi 1^ = ^Brcoo - Px (^'s„co being the Brcco 

four-momentum), and P+. can be calculated. The loss 
of one or more charged or neutral particles or the addi- 
tion of tracks or single electrons from photon conversions 
degrade the reconstruction of X„ and the resolution of 
the measurement of any related kinematic variables. In 
order to reduce the impact of missing charged particles 
and the effect of single electrons from 7 — >■ e+e" conver- 
sions, we impose charge conservation on the whole event, 
Qtot = Q-Broco + Qx + Qe ~ 0- This requirement rejects 
a larger fraction of i? — > XciD events because of their 
higher charged multiplicity and the presence of very low 
momentum charged pions from D*^ — > D'^T^^^f^ decays 
which have low detection efiiciency. 

In i? — > X£D decays, where the state X decays hadron- 
ically, the only undetected particle is a neutrino. The 
neutrino four-momentum P^, can be estimated from the 
missing momentum four- vector Pmiss = Pr{4S) ~ ^Srcco ^ 
Px — Pt ■ For correctly reconstructed events with a single 
semileptonic decay, the missing mass squared, MM'^ = 
P^iss, is consistent with zero. Failure to detect one or 
more particles in the event creates a tail at large positive 
values; thus MM^ is used as a measure of the quality of 
the event reconstruction. Though MM^ is Lorentz in- 
variant, the missing momentum is usually measured in 
the laboratory frame, because this avoids the additional 
uncertainty related to the transformation into the cm. 
frame. We require MM^ to be less than 0.5 GeV^. Be- 
cause of the higher probability for additional unrecon- 
structed neutral particles, a neutrino or Kl, the MM^ 
distribution is broader for B — > Xciiy decays, and this 
restriction suppresses this background more than signal 
events. 

In addition, we suppress the B — > D*£i> background 
by exploiting the small Q-value of the D* — )■ DtTsoH de- 
cays which result in a very low momentum pion. For 



energetic D* mesons, the momenta PTr^ott ^^"^ Pd are 
almost collinear, and we can approximate the D* di- 
rection by the TTgoft direction and estimate the D* en- 
ergy by a simple approximation based on the iJ-n-aoft i 
« mu* xi5T^^fj/145 MeV. Using the measured iJrcco 
and charged lepton momenta, and the four-momentum of 
the D* derived from any pion with cm. momentum be- 
low 200 MeV, we estimate the neutrino mass for a poten- 
tial B^ D*ep decay as MM^^^ {Pb-Pd* -Ptf- For 
true B — > D*iv decays, this distribution peaks at zero. 
Thus, we veto D* decays to low momentum charged or 
neutral pions by requiring, respectively, MM'^^^^{-K'^^f^) < 
-3 GeV2 or MM^^t„(7r0^j.J < -2 GeV^. This is achieved 
without explicit reconstruction of the D meson decays, 
and thus avoids large losses in rejection power for this 
veto. 

We reduce B — > D*£i> background by vetoing events 
with a charged or neutral kaon (K^ — >■ tt+tt^), that orig- 
inate primarily from the decays of charm particles. 

A summary of the impact of the signal selection cri- 
teria on the high-energy lepton sample, for the signal, 
semileptonic and nonsemileptonic background samples 
is presented in Table II, in terms of cumulative selec- 
tion efficiencies. Figure 1 shows the kinematic vari- 
ables that appear in Table II for different event cate- 
gories. Combinatorial background is not included; it is 
subtracted based on fits to the tties distributions, as de- 
scribed in Section III D. The overall efficiency for select- 
ing charmless semileptonic decays in the sample of tagged 
events with a charged lepton is 33.8%; the background 
reduction is 97.8% for B XJP and 95.3% for "other". 

The resolution functions determined from MC simula- 
tion of signal events passing the selection requirements 
are shown in Fig. 2 for the variables Mx, q^, and P+. 
Each of these distributions has a narrow core contain- 
ing 30%, 50%, and 30% of the B XuiD events, with 
widths of 25 MeV, 250 MeV^, and 10 MeV, respectively. 
The remaining events have a considerably poorer resolu- 
tion, primarily because of lost secondary particles from 
the decay of the hadronic X^. 

TABLE II: Comparison of the cumulative selection efficien- 
cies for samples of signal B — >■ Xuii' decays and B — >■ Xcii' 
and "other" backgrounds. The efficiencies are relative to the 
sample of Srcco-tagged events with a charged lepton. 



Selection 


B X^lv B 


-5. Xdv Other 


Only one lepton 


99.3% 


98.1% 


95.8% 


Total charge Q=0 


65.5% 


52.9% 


49.1% 




44.2% 


17.8% 


17.8% 


D*£u (7r+) veto 


40.6% 


9.9% 


14.4% 


D*£u {tv°) veto 


34.8% 


6.3% 


9.1% 


Kaon veto 


33.8% 


2.2% 


4.7% 



On the basis of the kaon and the D* veto, two data 
samples are defined: 

• signal- enriched: events that pass the vetoes; this 
sample is used to extract the signal; 
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FIG. 1: MC distribution of the kinematic variables for wliich we apply restrictions sequentially as listed in Table 11, for 
B — ^ Xuly (solid line) , B — ^ Xctv (dashed line) and "other" component (dotted line) . All distributions are normalized to unity 
and selection criteria have been applied cumulatively, except those affecting directly the variable shown. The arrows indicate 
the selection requirement for a specific variable, as described in Section IIIC. 



• signal- depleted: events rejected by at least one veto; 
they are used as control sample to check the agree- 
ment between data and simulated backgrounds, in- 
cluding the poorly understood B — > D**iv decays. 



D. Subtraction of combinatorial background 

The subtraction of the combinatorial background of 
the i?reco tag for the signal and normalization samples 
relies on unbinned maximum-likelihood fits to the mES 
distributions. For signal decays the goal is to extract the 
distributions in the kinematic variables p^, Mx^ 9^, and 
P+. Because the shapes and relative yields of the signal 



and background contributions depend on the values of 
these kinematic variables, the continuum and combina- 
torial background subtraction is performed separately for 
subsamples corresponding to events in bins of these vari- 
ables. This results in more accurate spectra than a single 
fit to the full sample of events in each selected region of 
phase space. 

For the normalization sample, the fit is performed for 
the full event sample, separately for and B^ tags. 

The Toes distribution for the combinatorial i?roco back- 
ground can be described by an ARGUS function [52] , 

/bkg(TO) = 7VbkgTOVl-m2e-«(i-™'), (3) 
where to = toes/toes™''^ and toes™'*'^ is the endpoint 
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FIG. 2: Resolution for MC simulated for signal B Xuii^ events passing all event selection criteria, (left) Mxr 



(center) 



and (right) P+ 



The curve shows a fit results for the sum of two Gaussian functions. 



of the TTiEs distribution which depends on the beam en- 
ergy, and ^ determines the shape of the function. A^bkg 
refers to the total number of background events in the 
distribution. 

For signal events, the mss distribution resembles a 
resolution function peaking at the B meson mass with 
a slight tail to lower masses. Usually the peak of the 
rriES distribution is empirically described by a Crystal 
Ball function [53], but this ansatz turned out to be in- 
adequate for this dataset because the Breco sample is 
composed of many individual decay modes with different 
resolutions. We therefore follow an approach previously 
used in BaBAR data [54] and build a more general func- 
tion, using a Gaussian function, fg{x) = and the 
derivative of tanhx, ft{x) = e~^/(l -I- e~^), to arrive at 



(C3-A)" 
Ci f / A N 



l-r 

0"2 



if A < a 

if a A < (4) 
if A ^ 0. 



Here A = tties ~ ^es, where mss is the maximum of 
the 771ES distribution. C\ , C2 and C3 are functions of the 
parameters mEs, ci, (T2, cl, a, and n, that ensure the 
continuity of Uig- 

Given the very large number of parameters, we first 
perform a fit to samples covering the full kinematic range 
and determine all parameters describing /^^g and the AR- 
GUS function. We then repeat the fit for events in each 
bin of the kinematic variables, with only the relative nor- 
malization of the signal and background, and the shape 
parameter ^ of the ARGUS function as free parameters. 
Figure 3 shows the toes distribution for the inclusive 
semileptonic sample, separately for charged and neutral 
B mesons. 

Finally, we correct for the contamination from cascade 
background in the number of neutral B mesons, due to 
the effect of B'^-B^ mixing, in each bin of the kinematic 
variables. We distinguish neutral B decays with right- 
and wrong-sign leptons, based on the flavor of the i?rcco 
decay. The contribution from cascade decays is sub- 
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FIG. 3: The tties distribution for the inclusive semileptonic 
sample, for fully reconstructed hadronic decays of _B~ (left) 
and B*^ mesons (right). The solid line shows the result of 
the maximum-likelihood fit to signal and combinatorial back- 
grounds; the dashed line indicates the shape of the back- 
ground described by an ARGUS function. 



tracted by computing the number of neutral B mesons 
Nbo as 



Xd 



2xd 1 - 2xrf 

where Nbo and Nbo are the number of neutral B mesons 
with right and wrong sign of the charge of the accompa- 
nying Icpton, and Xd = 0.188 ± 0.002 [40] is the B°-B" 
mixing parameter. 

The performance of the toes fit has been verified using 
MC simulated distributions. We split the full sample in 
two parts. One part, containing one third of the events, 
is treated as data, and is similar in size to the total data 
sample. The remaining two thirds represent the simu- 
lation. The fit procedure, described in Section IV, is 
applied to these samples and yields, within uncertain- 
ties, the charmless semileptonic branching fraction that 
is input to the MC generation. 
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IV. SIGNAL EXTRACTION AND PARTIAL 
BRANCHING FRACTION MEASUREMENT 

A. Signal yield 

Once continuum and combinatorial BB backgrounds 
have been subtracted and the mixing correction has been 
appUed, the resulting differential distributions of the 
kinematic variables are fitted using a minimization 
to extract Nu, the number of selected signal events. The 
for these fits is defined as 



(6) 



where, for each bin i of variable width, is the num- 
ber of observed events, and N^;^'-' and N^^'~^ are the 
number of MC predicted events for signal and back- 
ground, respectively. The statistical uncertainties (j{N^) 
and (j{N^'^^) are are taken from fits to the toes distribu- 
tions in data and MC simulations. The scale factors Csig 
and Cbkg are free parameters of the fit. The differential 
distributions arc compared with the sum of the signal and 
background distributions resulting from the fit in Figs. 4 
and 5. For the B — >■ Xu^v signal contributions we distin- 
guish between decays that were generated with values of 
the kinematic variable inside the restricted phase space 
regions, and a small number of events, -/V°"*, with values 
outside these regions. This distinction allows us to re- 
late the fitted signal yields to the theoretical calculations 
applied to extract \Vub\- 



B. Partial branching fractions 

We obtain partial branching fractions for charmless 
semilcptonic decays from the observed number of signal 
events in the kinematic regions considered, after correc- 
tion for background and efficiency, and normalization to 
the total number of semileptonic decays B — )■ Xtv ob- 
served in the -Broco event sample. For each of the re- 
stricted regions of phase space under study, we calculate 
the ratio 



Ai?u/sl 



Ag(B Xulv) _ 
B{B Xii>) 

{N^i-BGsi) e^L 



(7) 



of phase space, and ej* refers to the efficiency for select- 
ing a lepton from a i? — >■ Xuiv decay with a momentum 
P(>1 GeV in a signal event tagged with efficiency e"ag- 
Similarly, iV*™'' is related to Ns\, the fitted number of 
observed Breco accompanied by a charged lepton with 
p* > 1 GeV, through iV,Y"° = i^si- BGsi)/efelig. Here, 
BGsi is the remaining peaking background estimated 
from MC simulation and A^si is obtained from the toes 
fit to the selected semileptonic sample and e^' refers to 
the efficiency for selecting a lepton from a semileptonic B 
decay with a momentum p| > 1 GeV in an event tagged 
with efficiency ell^. We obtain N^i = 237, 433 ± 838 and 
BGsi = 20, 705 ± 132. 

The ratio of efficiencies in Eq. (7) accounts for dif- 
ferences in the final states and the different lepton mo- 
mentum spectra for the two classes of events, and their 
impact on the tagging. The efficiencies for Srcco tagging 
and lepton detection are not very different, and thus the 
efficiency ratio is close to one. 

We convert Eq. (7) to partial branching fractions by 
using the total semileptonic branching fraction, B{B — > 
X£P) = (10.75 ±0.15)% [40]. 

The regions of phase space, fitted event yields, efficien- 
cies introduced in Eq. (7), and partial branching fractions 
are listed in Table III; the regions are one-dimensional in 
Mx, P+-, or p^, or two-dimensional in the plane Mx ver- 
sus q^. In the following, we will refer to the latter as 
Mx - ■ Two fits have been performed with no addi- 
tional kinematic restrictions, apart from the requirement 
p*f > \ GeV: a fit to the lepton momentum spectrum and 
a fit to the two-dimensional histogram Mx - ■ Since 
the same events enter both fits, the correlation is very 
high. The fact that the results arc in excellent agreement 
indicates that the distribution of the simulated signal and 
background distributions agree well with the data. 

Correlations between the different analyses are re- 
ported in the entries above the main diagonal of Table IV. 

In addition, a series of fits to the lepton momentum 
spectrum has been performed with the lower limit on 
p\ increasing from 1.0 GeV to 2.4 GeV. The results 
are presented in Section VI; the measurement at p*^ > 
1.3 GeV gives the smallest total uncertainty and is also 
quoted in Table III. 

Consistency checks have been performed. The analysis 
done on data samples collected in different data-taking 
periods, or separating the lepton flavor or charge have all 
yielded the same results, within experimental uncertain- 
ties. 



Here, iV^''"° and N^™° refer to the true number of signal 
and normalization events. The observed signal yield 
is related to iV*--- by iV„ = el,eJ,^e]tel^,^N:,^-^ , where e^, 
is the efficiency for detecting B — >■ XuiP decays in the 
tagged sample after applying all selection criteria, ekin is 
the fraction of signal events with both true and recon- 
structed Mx, -P+, 9^, or within the restricted region 



C. Partial branching fractions for and B 

All the fits, except those to the p| distribution, have 
been repeated separately for charged and neutral Breco 
tags. In this case, we extract the true signal yields from 
the measurements by the following relations to determine 
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the partial branching fractions: 



where the cross-feeds probabilities, V^- _^ga and 
'Pb" -i-B" ' ^'"^ computed using MC simulated events 
and are typically of the order of (2 - 3)%. 

Figure 6 shows the distributions of B X.y_£iy 
events after background subtraction, for charged and 
neutral B decays, with Mx < 1.7 GeV. Fitted yields, 
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TABLE III: List of the fitted numbers of signal events Nu, the number of events generated outside the kinematic selection 
N°^^ , the efficiencies, the partial branching fractions AB{B Xutv) and the per degree of freedom for the different selected 
regions of phase space. The first uncertainty is statistical, the second systematic. The > 1 GeV requirement is implicitly 
assumed. 



Region of phase space 




»rout 


u u 
scl kin 




AB{B Xulv) {IQ-'-') 


xVndof 


Mx < 1.55 GeV 


1033 ± 73 


29 ± 2 


0.365 ± 0.002 


1.29 ±0.03 


1.08 ±0.08 ±0.06 


7.9/8 


Mx < 1.70 GeV 


1089 ± 82 


25 ± 2 


0.370 ± 0.002 


1.27 ±0.04 


1.15 ±0.10 ±0.08 


6.6/8 


P+ < 0.66 GeV 


902 ± 80 


54 ±5 


0.375 ± 0.003 


1.22 ±0.03 


0.98 ±0.09 ±0.08 


3.4/9 


Mx < 1.70 GeV, > 5 


! GeV^ 665 ±53 


39 ±3 


0.386 ± 0.003 


1.25 ±0.03 


0.68 ±0.06 ±0.04 


23.7/26 


Mx-q^ 


1441 ± 102 





0.338 ± 0.002 


1.18 ±0.03 


1.80 ±0.13 ±0.15 


31.0/29 


Pt > 1.0 GeV 


1470 ± 130 


8±2 


0.342 ± 0.002 


1.18 ±0.03 


1.81 ±0.16 ±0.19 


21.6/14 


Pe > 1.3 GeV 


1329 ± 121 


61 ±5 


0.363 ± 0.002 


1.18 ±0.09 


1.53 ±0.13 ±0.14 


20.4/14 



TABLE IV: Correlation coefficients for measurements in different kinematic regions. The entries above the main diagonal refer 
to correlations (statistical and systematic) for pairs of measurements of the partial branching fractions; the entries below the 
diagonal refer to the correlations (experimental and theoretical) for pairs of \Vub\ measurements. 



Phase space restriction 


Mx < 1.55 
GeV 


Mx < 1.70 
GeV 


P+ < 0.66 
GeV 


Mx < 1.70 GeV, 
> 8 GeV^ 


Mx -q" 
p'e > 1.0 GeV 


p1 > 1.3 
GeV 


Mx < 1.55 GeV 


1 


0.77 


0.74 


0.50 


0.72 


0.57 


Mx < 1.70 GeV 


0.81 


1 


0.86 


0.55 


0.94 


0.73 


P+ < 0.66 GeV 


0.69 


0.81 


1 


0.46 


0.78 


0.61 


Mx < 1.70 GeV, > 8 GeV^ 


0.40 


0.46 


0.38 


1 


0.52 


0.46 


Mx~q^ 


0.58 


0.88 


0.67 


0.34 


1 


0.74 


p| > 1.3 GeV 


0.53 


0.72 


0.58 


0.40 


0.72 


1 



efBciencies, and partial branching fractions are given in 
Table V. 



D. Data - Monte Carlo comparisons 

The separation of the signal events from the noncom- 
binatorial backgrounds relies heavily on the MC simula- 
tion to correctly describe the distribution for signal and 
background sources. Therefore, an extensive study has 
been devoted to detailed comparisons of data and MC 
distributions. 

A correction applied to the simulation improves the 
quality of the fits to the kinematic distributions in re- 
gions that are dominated hy B Xc(-v background, es- 
pecially in the high Mx region. In the simulation, we 
adjust \d", the ratio of branching fractions of semilep- 
tonic decays to P-wave D mesons and nonresonant charm 
states decaying to D^*^X, over the sum of all D'^*hv and 
"other" background components, 

^ _ BiB ^ D**iD) +B{B ^ D^*'>X£D) 
^ B{B~^Di*)lD)+B{B^X,t^,,) ■ 

This ratio has been determined from data by performing 
a fit on the Mx - distribution of the signal-depleted 
sample without kinematic selection. The resulting dis- 
tribution of this fit is shown in Fig. 7. We measure 



Ad" = 0.73 ± 0.08, where the error takes into account 
the fact that x^/ndof = 2. Other determinations, using 
signal-enriched samples, give statistically consistent re- 
sults. This adjustment improves the quality of the fits in 
regions where backgrounds dominate, but it has a small 
impact on the fitted signal yield. We have verified that 
using D** MC correction factors determined separately 
on each analysis do not change significantly the results 
with respect to our default strategy, where A^*. is deter- 
mined for the most inclusive sample available, namely the 
signal-depleted sample of the analysis without kinematic 
requirements. 

Figures 8 and 9 show comparisons of data and MC dis- 
tributions, after subtraction of the combinatorial back- 
ground, for signal-enriched and signal-depleted event 
samples. All the selection criteria have been applied, 
except those affecting directly the variable shown. The 
spectra are background-subtracted based on the results 
of the TTT-Es fit performed for each bin of the variable 
shown. The uncertainties on data points are on the yields 
of the bin-by-bin fits. The data and MC distributions are 
normalized to the same area. The overall agreement is 
reasonable, taking into account that the uncertainties are 
purely statistical. The effects that introduce differences 
between data and simulation are described in Section V; 
their impact is assessed and accounted for as systematic 
uncertainty. 
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TABLE V: Summary of the fits to separate samples of neutral and charged B decays. For details see Table III. 



B" decays 








u u 
^sol ^kin 




AB{B X^lu) (10-3) 


X^/ndof 


Mx < 1.55 GeV 




458 ± 48 


12 ± 1 


0.360 ± 0.004 


1.49 ± 0.07 


1.09 ± 0.12 ± 0.11 


19.0/9 


Mx < 1.70 GeV 




444 ± 53 


12 ± 1 


0.370 ± 0.004 


1.45 ± 0.07 


1.12 ± 0.11 ± 0.11 


16.6/9 


P+ < 0.66 GeV 




434 ± 52 


27 ± 3 


0.367 ± 0.004 


1.38 ± 0.06 


1.09 ±0.13 ±0.11 


9.1/9 


Mx < 1.70 GeV, > 5 


GeV^ 


262 ± 38 


16 ±2 


0.380 ± 0.005 


1.43 ±0.06 


0.61 ±0.09 ±0.06 


15.8/26 


Mx~q^ 




553 ± 72 





0.328 ± 0.003 


1.36 ±0.08 


1.58 ±0.21 ±0.20 


14.8/29 


B~ decays 






AfOUt 

u 


■u u 
scl kin 




AB(B^X^lu) (10-3) 


xVndof 


Mx < 1.55 GeV 




591 ± 56 


17 ±2 


0.370 ± 0.003 


1.18 ±0.04 


1.12±0.11±0.11 


3.1/9 


Mx < 1.70 GeV 




669 ± 63 


14 ± 1 


0.370 ± 0.003 


1.17 ±0.07 


1.27 ±0.14 ±0.13 


3.3/9 


P+ < 0.66 GeV 




491 ±61 


28 ±4 


0.379 ± 0.004 


1.11 ±0.03 


0.96 ±0.12 ±0.12 


2.0/9 


Mx < 1.70 GeV, > 5 


GeV^ 


406 ± 41 


24 ±2 


0.392 ± 0.004 


1.43 ±0.03 


0.74 ± 0.08 ± 0.08 


26.9/26 


Mx~q^ 




859 ± 79 





0.345 ± 0.003 


1.07 ±0.03 


1.91 ±0.18 ±0.22 


36.7/29 



CM 140 




FIG. 6: Comparison of the measured q^ distributions (data 
points) for Mx < 1.7 GeV for charmless semileptonic decays 
of (a) charged and (b) neutral B mesons to the results of 
the fit (histogram), after B — ^ X^iv and "other" background 
subtraction. 

V. SYSTEMATIC UNCERTAINTIES 

The experimental technique described in this article, 
namely the measurement of a ratio of branching fractions, 
ensures that systematic uncertainties due, for example. 



to radiative corrections or differences between B and 
or B^ production rate and lifetime, are negligible. 
A summary of all other statistical and systematic un- 
certainties on the partial branching fractions for selected 
kinematic regions of phase space is shown in Table VI for 
the complete data sample, and in Table VII for charged 
and neutral B samples separately. 

The individual sources of systematic uncertainties are, 
to a good approximation, uncorrelated and can therefore 
be added in quadrature to obtain the total systematic 
uncertainties for partial branching fraction. In the fol- 
lowing, we discuss the assessment of the systematic un- 
certainties in detail. 

To estimate the systematic uncertainties on the ratio 
AR^/gi , we compare the results obtained from the nom- 
inal fits with results obtained after changes to the MC 
simulation that reflect the uncertainty in the parameters 
which impact the detector efficiency and resolution or 
the simulation of signal and background processes. For 
instance, we lower the tracking efficiency by randomly 
eliminating a fraction of tracks (corresponding to the es- 
timated uncertainty) in the MC sample, redo the event 
reconstruction and selection on the recoil side, perform 
the fit, and take the difference compared to the results 
obtained with the nominal MC simulation as an estimate 
of the systematic uncertainty. The sources of systematic 
uncertainties are largely identical for all selected signal 
samples, but the size of their impact varies slightly. 



A. Detector effects 

Uncertainties in the reconstruction efhciencies for 
charged and neutral particles, in the rate of tracks 
and photons from beam background, misreconstructed 
tracks, failures in the matching of EMC clusters to 
charged tracks, showers split-off from hadronic interac- 
tions, undetected i^^, and additional neutrinos, all con- 
tribute to the event reconstruction and impact the vari- 
ables that are used in the event selection and the analysis. 
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FIG. 7: Fit results to the Mx - distribution for the signal-depleted sample. The distribution is reported separately for 
the four Mx bins: (a) Mx < 1.5 GeV, (b)_1.5 < Mx < 2.0 GeV, (c) 2.0 < Mx < 2.5 GeV and (d) 2.5 < Mx <_3.0 GeV. The 
three MC contributions shown here are: B — X„£u decays vetoed by the selection (no shading), B — s- D£u, B — > D*£u and 
"other" background (light shading), and the B — >■ D**£ii component as defined in the text (dark shading). 



For all these effects the uncertainties in the efficiencies 
and resolution have been derived from comparisons of 
data and MC simulation for selected control samples. 

From the study of the angular and momentum distri- 
butions of low momentum pions in D* samples, we es- 
timate the uncertainty on the track finding efficiency at 
low momenta to be about 1.0%. For all other tracks, the 
difference between data and MC in tracking efficiency is 
estimated to be about 0.5% per track. The systematic 
uncertainty on the ratio Ai?„/^; is calculated as described 
above, and shown in Tables VI and VII. 

Similarly, for single photons, we estimate the system- 
atic uncertainty by randomly eliminating showers that 
are not matched to the Tr^^fj used to veto B D*iD 
decays, with a probability of 1.8% per shower. 

We estimate the systematic uncertainty due to 7r° de- 
tection by randomly eliminating neutral pions that are 
used in the B — )- D*tv veto, with a probability of 3% per 

Uncertainties on charged particle identification effi- 
ciencies have been assessed to be 2.0% for electrons and 
3.0% for muons. The uncertainty on the corresponding 
misidentification rates are estimated to be 15%. System- 
atic uncertainties on the kaon identification efficiency and 
misidentification rate are 2% and 15%, respectively. 

In this analysis, no effort was made to identify K^. On 
the other hand, mesons interacting in the detector 
deposit only a fraction of their energy in the EMC, thus 
they impact Pmiss and other kinematic variables used in 



this analysis. Based on detailed studies of data control 
samples of — >■ if^Tr+Tr" decays, corrections to the if^ 
efficiency and energy deposition have been derived and 
applied to the simulation as a function of the K'^ momen- 
tum and angles. We take the difference compared to the 
results obtained without this correction applied to the 
simulation as an estimate of the systematic uncertainty. 

Differences in both and A'^ production rates of data 
and MC are taken into account by adjusting the inclusive 
D -> K^X and Ds -> K^X branching fractions. The 
associated systematic uncertainty is assessed by varying 
these branching fractions within their uncertainties. 



B. Signal and background simulation 

1. Signal simulation 

Knowledge of the details of inclusive B X^tv decays 
is crucial to several aspects of the analysis: the fraction 
of events within the selected kinematic region depends on 
the signal kinematics over the full phase space. Specif- 
ically, the efficiencies and ekin I'dy on accurate MC 
simulation, because the particle multiplicities, momenta, 
and angles depend on the hadronization model for the 
hadronic states X„. 

To simulate the signal B — > X^tv decays we have cho- 
sen the prescription by De Fazio and Neubert [38] . Differ- 
ent choices of the parameterization for the Fermi motion 
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FIG. 8: Comparison of data (points with statistical uncertainties) and MC (histograms) simulated distributions of (a,d) the 
missing mass squared, (b,e) the missing momentum, and (c,f) the missing energy for B — >■ Xuiv enhanced (top row) and 
depleted (bottom row) event samples. 



of the b quark inside the B meson (Section II D) lead to 
different spectra of the hadron mass Mx and lepton mo- 
mentum . We estimate the impact of these choices by 
repeating the analysis with shape function parameters set 
to values of Af^and A^^ corresponding to the contour 
of the Ax^ — 1 error ellipse [55]. To assess the impact of 
the choice of the SF ansatz, we repeat this procedure for 
a different SF ansatz [38]. 

Since the simulation of i? — )• Xutv decays is a hybrid 
of exclusive decays to low-mass charmless mesons and 
inclusive decays to higher-mass states Xu, the relative 
contributions of the various decays impact the overall 
kinematics and thereby the efficiencies. Wc evalute the 
impact of varying the branching fractions of the exclu- 
sive charmless semileptonic B decays by one standard 
deviation. 

The signal losses caused by the kaon veto depend on 
the production rate of kaons in these decays. In the MC 
simulation, the number of A'^ and A'^ in the signal de- 
cays is set by the probability of producing ss quark pairs 
from the vacuum. The fraction of ss events is about 
12.0% for the nonresonant component of the signal and 
is fixed by the parameter 7., in the fragmentation by 
JETSET [39]. This parameter has been measured by two 



experiments at center of mass energies between 12 and 
36 GeV as 7^ = 0.35 ± 0.05 [56], 7^ = 0.27 ± 0.06 [57]. 
We adopt the value 7s = 0.3 and estimate the systematic 
uncertainty by varying the fraction of ss events by ±30%. 

The theoretical uncertainty due to the lower limit 
on the lepton spectrum is largely accounted for by the 
rcweighting of events for the assessment of the theoreti- 
cal uncertainty related to the Fermi motion. 



2. Branching fractions for B and D decays 

The exclusive semileptonic branching fractions for 
B Xclv decays and the hadronic mass spectra for 
these decays are crucial for the determination of the yield 
of the inclusive normalization sample and the B — > Xciv 
background. Exclusive B and D branching fractions used 
in the MC simulation differ slightly from the world av- 
erages [40]; this difference is corrected by reweigliting 
events in the simulation. The branching fraction for 
the sum of semileptonic decays to nonresonant D^*^tt 
or broad D** states is taken as the difference between 
the total semileptonic rate and the other well measured 
branching fractions, and amounts to about 1.7%. 
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FIG. 9: Comparison of data (points with statistical uncertainties) and MC (histograms) simulated distributions of (a,d) the 
charged track multiplicity, (b,e) the photon multiplicity, (c,f ) and the total charge per event for B — >■ Xulv enhanced (top row) 
and depleted (bottom row) event samples. 



Similarly, branching fractions and decay distributions 
for hadronic and semileptonic D meson decays affect the 
measurement of Ai?^/^; . The effect is different for neutral 
and charged B mesons, because decays mostly into 
charged D mesons while B~ decays almost always into 
neutral charm mesons. 

Likewise, uncertainties on the form factors for B — 
D^*^ii' decays are taken into account by repeating the 
analysis with changes of the form factor values by their 
experimental uncertainties [47]. For B — >■ D**tv decays, 
the uncertainties on the form factor have not been spec- 
ified. Thus, we perform the fits with the ISGW2 [58] pa- 
rameterization of the form factors and take the difference 
with respect to the default fits as systematic uncertainty. 

The uncertainty related to the Ad** parameter intro- 
duced in Eq. (8) has been estimated by varying it within 
its uncertainty, and taking the difference with respect to 
the default fits as systematic uncertainty. 

3. Combinatorial background subtraction and normalization 

For the fits to the toes distributions in individual bins 
of a given kinematic variable, all parameters other than 



event yields and the ARGUS shape are fixed to values de- 
termined from distributions obtained from the full signal 
sample. To estimate the systematic uncertainty due this 
choice of parameters, their values are varied within their 
statistical uncertainty, taking correlations into account. 
We estimate the effect of the combinatorial background 
subtraction by determining it on a simulated sample by 
means of Monte Carlo truth information, and getting the 
signal yields on data by subtraction. The differences rel- 
ative to the default fit are taken as systematic uncertain- 
ties. 

Finally, the uncertainty on the knowledge of the total 
semileptonic branching fraction adds 1.4% to the assess- 
ment of our systematic uncertainty. 

In summary, the smallest statistical and systematic un- 
certainties are achieved for the Mx < 1.55 GeV region, 
which has an acceptance that is reduced by 40% with re- 
spect to the region defined by > 1.0 GcV, but has the 
best separation of signal and background. The dominant 
systematic uncertainty for samples with no phase space 
restrictions, except for > 1.0 GcV, is due to the un- 
certainty on the shape function parameters which impact 
the differential and distributions. 
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TABLE VI: Statistical and systematic uncertainties (in percent) on measurements of the partial branching fraction in seven 
selected kinematic regions. The total systematic uncertainty is the sum in quadrature of the MC statistical uncertainty and all 
the other single contibutions from detector effects, signal and background simulation, background subtraction and normalization. 
The total uncertainty is the sum in quadrature of the data statistical and total systematic uncertainties. 



Phase space restriction 


Mx < 1.55 
GeV 


Mx < 1.70 P+ < 0.66 Mx 
GeV GeV 


< 1.70 GeV, 
> 8 GeV^ 


Mx - q^ 


Pi > 1.0 
GeV 


pI > 1.3 
GeV 


Data statistical uncertainty 


7.1 


8.9 


8.9 


8.0 


7.1 


9.4 


8.8 


MC statistical uncertainty 


1.3 


1.3 


1.3 


1.6 


1.1 


1.1 


1.2 


Detector effects 


Track efficiency 


0.4 


1.0 


1.1 


1.7 


0.7 


1.2 


1.0 


Photon efficiency 


1.3 


2.1 


4.0 


0.7 


1.0 


0.9 


0.9 


tt" efficiency 


1.2 


0.9 


1.1 


0.9 


0.9 


2.9 


1.1 


Particle identification 


1.9 


2.4 


3.3 


2.9 


2.3 


2.9 


2.2 


Kl production/detection 


0.9 


1.3 


1.1 


2.1 


1.6 


1.3 


0.6 


Ks production/detection 


0.8 


1.4 


1.7 


2.1 


1.2 


1.3 


0.3 






Sif 


;nal simulation 










Shape function parameters 


2.0 


1.3 


1.2 


0.7 


5.4 


6.4 


6.6 


Shape function form 


1.2 


1.6 


2.6 


1.2 


1.5 


1.1 


1.1 


Exclusive B Xuii' 


0.6 


1.3 


1.6 


0.7 


1.9 


5.3 


3.4 


ss production 


1.2 


1.6 


1.1 


1.0 


2.7 


3.1 


2.4 






Backj 


ground simulation 










B semileptonic branching ratio 


0.9 


1.4 


1.5 


1.4 


1.0 


0.8 


0.7 


D decays 


1.1 


0.6 


1.1 


0.6 


1.1 


1.6 


1.5 


B — )> Dlu form factor 


0.5 


0.5 


1.3 


0.4 


0.4 


0.1 


0.2 


B — >■ D*lv form factor 


0.7 


0.7 


0.9 


0.7 


0.7 


0.7 


0.7 


B — > D**lu form factor 


0.8 


0.9 


1.3 


0.4 


0.9 


1.0 


0.3 


B — ^ D** reweighting 


0.5 


1.4 


1.5 


1.0 


1.9 


0.4 


1.5 


rriES background subtraction 


m_BS background subtraction 


2.0 


2.7 


1.9 


2.6 


1.9 


2.0 


2.5 


combinatorial backg. 


1.8 


1.8 


2.6 


1.8 


1.0 


2.1 


0.5 


Normalization 


Total semileptonic BP 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


Total systematic uncertainty 


5.5 


6.7 


8.3 


6.6 


8.4 


11.0 


9.3 


Total experimental uncertainty 


9.0 


11.1 


12.2 


10.4 


11.0 


14.4 


12.8 



VI. EXTRACTION OF \Vnb\ 
A. QCD corrections 

We extract \Vub \ from the measurements of the partial 
branching fractions AB{B Xuiv) by relying on QCD 
predictions. In principle, the total rate for B — >■ Xuiv 
decays can be calculated based on heavy quark expan- 
sions (HQE) in powers of 1 /nib with uncertainties at the 
level of 5%, in a similar way as for B — >■ Xclv decays. 
Unfortunately, the restrictions imposed on the phase 
space to reduce the large background from Cabibbo- 
favored decays spoil the HQE convergence. Perturba- 
tive and nonperturbative corrections are drastically en- 
hanced and the rate becomes sensitive to the Fermi 
motion of the b quark inside the B meson, introduc- 
ing terms that are not suppressed by powers of l/rrib. 
In practice, nonperturbative SFs are introduced. The 
form of the SFs cannot be calculated from first princi- 
ples. Thus, knowledge of these SFs relies on global fits 
performed by several collaborations to moments of the 



lepton energy and hadronic invariant mass in semilep- 
tonic B decays, and of the photon energy in radiative 
B Xg^ inclusive decays [59-61]. We adopt results of 
the global fits to published measurements of moments, 
performed in the kinetic renormalization scheme, specif- 
ically the h quark mass m^"' = (4.560 ± 0.023) GeV 
and the mean value of the b quark momentum operator 
^2(Mn) ^ (0.453 ±0.036) GeV2 [46, 55]. Due to confine- 
ment and nonperturbative effects the quantitative values 
of the quark mass and other HQE parameters are specific 
to the theoretical framework in which it is defined. Thus 
the results of the global fits need to be translated to other 
schemes, depending on the QCD calculation used to ex- 
tract \Vub\- In the following, we determine \Vub\ based on 
four different QCD calculations. The numerical calcula- 
tions are based on computer code kindly provided by the 
authors. 

The measured partial branching fractions Ai3(i? — > 
Xu^v) are related to \Vub\ via the following equation, 

y Tb theory 
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TABLE VII: Statistical and systematic uncertainties (in percent) on the partial branching fraction for neutral and charged B 
mesons for the five selected kinematic regions. The total systematic uncertainty is the sum in quadrature of the MC statistical 
uncertainty and all the other single contibutions from detector effects, signal and background simulation, background subtraction 
and normalization. The total uncertainty is the sum in quadrature of the the data statistical and total systematic uncertainties. 



Phase space restriction 


ivix 


1.00 

GeV 


IVIx 


<^ i. / U 

GeV 


p 


<. U.DD 

GeV 


Mx < 1.70 GeV, 
(7^ > 8 GeV^ 


Mx 






B° 


B' 


w 


B- 


B« 


B- 


B° 


B- 


w 


B- 


Data statistical uncertainty 


10.4 


9.6 


14.4 


11.0 


12.0 


12.5 


14.6 


10.1 


13.0 


9.2 


MC statistical uncertainty 


2.5 


1.6 


2.5 


1.6 


2.4 


1.8 


2.8 


2.0 


1.9 


1.3 


Detector effects 


4.5 


4.9 


5.0 


6.3 


5.9 


7.2 


5.3 


6.3 


5.6 


4.7 


Signal simulation 


6.6 


5.4 


5.2 


4.8 


4.8 


5.9 


3.7 


5.4 


8.7 


7.6 


Background simulation 


4.4 


4.2 


5.6 


4.5 


5.9 


5.4 


4.9 


4.4 


4.4 


5.0 


mES background subtraction 


4.1 


5.4 


5.2 


5.0 


2.9 


5.3 


5.2 


5.1 


3.8 


4.1 


Total semileptonic BP 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


1.4 


Total systematic uncertainty 


10.4 


10.2 


10.9 


10.6 


10.5 


12.2 


10.1 


11.0 


12.1 


11.2 


Total experimental uncertainty 


14.7 


14.0 


18.1 


15.3 


15.9 


17.5 


17.8 


14.9 


17.8 


14.5 



where AFthoory, the theoretically predicted B — )■ Xu^v 
rate for the selected phase space region, is based on dif- 
ferent QCD calculations, and the B meson hfetime is 
TB = 1.582 ± 0.007 ps [46]. We adopt the uncertain- 
ties on AFthoory as assessed by the authors. It should be 
noted that the systematic uncertainty on the branching 
fraction that is related to the uncertainties on the SF 
parameterization are fully correlated to the theoretical 
uncertainties discussed here. 

The calculated decay rates AFthoory and the resulting 
\Vub\ values are shown for the various kinematic regions 
in Tables VIII and IX, separately for the four different 
QCD calculations. 



1. BLNP calculation 

The theoretical uncertainties [18-20] arise from the un- 
certainty on TTib, /xj and other nonperturbative correc- 
tions, the functional form of the leading and the sub- 
leading SFs, the variation of the matching scales, and 
the uncertainty on the estimated contribution from weak 
annihilation processes. The dominant contributions are 
due to the uncertainties on mi,, and ^J. These parame- 
ters need to be translated to the shape function renormal- 
ization scheme, for which mf^'^ = (4.588 ± 0.025) GeV 

and M^^^^^ = (0.189 ±0.051) GeV^. The stated errors in- 
clude the uncertainties due to higher order terms which 
are neglected in the translation from one scheme to an- 
other. 

A recent calculation at NNLO [62] indicates that the 
differences with respect to NLO calculations are rather 
large. They would increase the value of \Vub\ by about 
8%, suggesting that the current uncertainties are under- 
estimated. Similar effects might also be present for other 
QCD calculations, but estimates are not yet available. 



2. DGE calculation 

The theoretical uncertainties [21, 22] arise from the un- 
certainty on as, the uncertainty on to;,, and other non- 
perturbative corrections, for instance, the variation of the 
matching scales, and the uncertainty on the weak annihi- 
lation. The dominant error is the uncertainty on to^ for 
which the MS renormalization scheme is used. There- 
fore the results of the global fit had to be translated to 
the MS scheme, = (4.194 ± 0.043) GeV, where the 
uncertainty includes the uncertainty on the translation. 

3. GGOU calculation 

The theoretical uncertainties [25] in the determinations 
of the widths and \ Vab\ from the GGOU calculations arise 
from the uncertainty on a^, to^, and plus various 
nonperturbative corrections: the modeling of the tail 
and choice of the scale gj, the functional form of the 
distribution functions, and the uncertainty on the weak 
annihilation rate. The dominant error originates from the 
uncertainties on mt and /ij. Since GGOU calculations 
are based on the kinetic renormalization scheme, there is 
no need for translation. 

4- ADFR calculation 

The ADFR calculation [23, 24] relates AS(B ^ X^lv) 
to \Vub\ in a way that is different from the other three 
calculations discussed above. In the framework of ADFR, 
the partial branching ratio is expressed in terms of Rc/m 

ABiB^XJ,)^^i^-^W, (10) 

1 + ^c/u 

where W = AT{B XJD)/T{B XJP) is the frac- 
tion of the charmless branching fraction in a selected 



kinematic region and B{B Xlv) is the total semilep- 
tonic branching fraction. Rc/u is related to \Vub\ as 

Rc,u^^-^I{p)G{as,p). (11) 

The function I{p) accounts for the suppression of phase 
space due to rric and I{p) = 1 — 8p+12p^log(l/p) + 8p^ — 
p"*, with p = rn^/ml w 0.1. The factor G{as,p) contains 
corrections suppressed by powers of and powers of p, 

oo 

G(a„p) = l + ^G„(pX, (12) 

n=l 

with G„(0) = 0. The errors of the ADFR calculations 
arise from the uncertainty in as, |Vcb|, the quark masses 
nih and mc, and the uncertainty on B{B — XiP). The 
dominant uncertainty is due to the uncertainty on the 
mass rric- 



B. \Vub\ extraction 

We present the results for \Vub\ with statistical, sys- 
tematic and theoretical uncertainties in Table VIII. Val- 
ues of \Vub\ extracted from partial branching fractions 
for samples with the lower limit on the lepton momen- 
tum pI varying from 1.0 GeV to 2.4 GeV are tabulated 
in Table IX. The different values of \Vub\ are consis- 
tent within one standard deviation and equally consistent 
with the previous BaBAR measurements of | Vub \ on inclu- 
sive charmless semileptonic B decays [11, 15, 16] as well 
as a similar measurement by the Belle Collaboration [17]. 

Our result on the study of the lepton spectrum above 
2 GeV can be compared to what BaBAR [15, 16], Belle [14] 
and CLEO [13] have published on the analysis of the lep- 
ton endpoint spectrum in untagged B decays. Experi- 
mental uncertainties are comparable, as well as theoret- 
ical uncertainties, which are quite large in this region of 
phase space. The values of \Vub\ obtained with such dif- 
ferent techniques agree very well. 

We have evaluated the correlations of the measure- 
ments of \Vub\ in selected regions of phase space tak- 
ing into account the experimental and theoretical pro- 
cedures, as presented in Table IV. The theoretical cor- 
relations have been obtained for the BLNP calculations 
by taking several values of the heavy quark parameters 
within their uncertainties and computing the correlation 
of the acceptance for pairs of phase space regions. The 
resulting correlation coefficients are in all cases greater 
than 97%. It is assumed that the correlations are also 
close to 100% for the other three theory calculations. 

We choose to quote the \ Vub\ value corresponding to the 
most inclusive measurement, namely the one based on 
the two-dimensional fit of the Mx - distribution with 
no phase space restrictions, except for p| > 1.0 GeV. 
We calculate the arithmetic average of the values and 
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uncertainties obtained with the different theoretical cal- 
culations shown above and find 

\Vub\ = (4.33 ± 0.24 ± 0.15) x 10"^ (13) 

where the first uncertainty is experimental and the sec- 
ond theoretical. 

A calculation specifically suited for phase space regions 
defined by the Mx and cuts [63] can also be consid- 
ered. This uses as input the b quark mass in the IS 
scheme [64], mj^ = 4.704 ± 0.029 GeV, determined by a 
global fit in that scheme, similar to the one described in 
Section VIA. The resulting value of \Vab\ for the phase 
space region defined by Mx < 1.7 GeV, > 8 GeV^ 
is \Vub\ = (4.50 ± 0.24 ± 0.29) x 10"'\ slightly larger but 
still in agreement with the other theoretical calculations. 

C. Limits on weak annihilation 

The measurements of AB{B — s- X^tv), separately for 
neutral and charged B mesons, are summarized in Ta- 
ble V for the various kinematic selections. These results 
are used to test isospin invariance, based on the ratio 

AT^ ^ AB{B- ^ XJv) 
Arc r" AB{B" ^ Xuif)' 

where r^/r" = 1.071 ± 0.009 [40] is the ratio of the life- 
times for B~ and B^. For the Mx < 1-55 GeV selection, 
we obtain R - 1 = 0.03 ± 0.15 ± 0.18, where the first 
uncertainty is statistical and the second is systematic. 
This result is consistent with zero; similar results, with 
larger uncertainties, are obtained for the other regions of 
phase space listed in Table V. Thus, we have no evidence 
for a difference between partial decay rates for B~ and 
If we define the possible contribution of the weak 
annihilation as ATwA ~ AF^ — AF'^, its relative con- 
tribution to the partial decay width AF for B — >■ Xu^P 
decays is ATwa/AT — R — 1. With fwA defined as 
the fraction of weak annihilation contribution for a spe- 
cific kinematic region and /„ defined as the fraction of 
B — > Xuii' events predicted for that region, we can write 
ATwA = fwA^WA and AF = /„F, where F is the to- 
tal decay width oi B ^ Xu^v decays. Thus the relative 
contribution of the weak annihilation is 

^^ = -^(i?-l). (15) 

1 JWA 

Since the weak annihilation is expected to be confined to 
the high q^ region, it is reasonable to assume fwA = 1.0 
for all the kinematic selections. We adopt the predic- 
tion for /„ by De Fazio-Neubert (see Section II D) and 
place limits on Twa/'^- The most stringent limit is 
obtained for the selection Mx < 1-55 GeV, namely 
-0.17 < (FvKA/r) < 0.19 at 90% confidence level (C.L.). 
This model-independent limit on WA is consistent, but 
weaker than the limit derived by the CLEO collabora- 
tion [65] on the basis of an assumed q^ distribution. 
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TABLE VIII: Results for \Vub\ obtained for the four different QCD calculations. The sources of the quoted uncertainties are 
experimental statistical, experimental systematic and theory, respectively. The theoretical B — >■ Xiu widths, AFthoory in ps^^, 
for the various phase space regions examined, as determined from the BLNP, DGE and GGOU calculations, are also shown. 
The ADFR calculation uses another methodology (see text), therefore the values for AFthcory have been obtained by inverting 
Eq. (9). The Pi > 1 GeV requirement is implicitly assumed in the definitions of phase space regions, unless otherwise noted. 



QCD Calculation Phase Space Region 



AF 



theory 



\Vub\{lO- 



BLNP 



Mx < 1.55 CeV 

Mx < 1.70 GeV 

P+ < 0.66 GeV 

Mx < 1.70 GeV, > 8 GeV^ 

Mx - g^ Pe > 1.0 GeV 

p't > 1.0 GeV 

p*f > 1.3 GeV 



H4.7 
-4.3 



39.3^ 
46.llf:o 
38.3+*" 



-4.3 
^3.0 
-2.4 
^6.2 
-5.0 

62.0+«-2 



23.1 
62.0! 



52. 



-5.0 
5+5.3 
^-4.3 



4.17 ±0.15 ±0.12^0 M 
3.97 ± 0.17 ±0.14to:2o 
4.02 ± 0.18 ±0.16l:;];23 
4.25 ± 0.19 ±0.13l:^;25 

4.28 ± 0.15 ±0.18t^;2o 
4.30 ± 0.18 ±0.2lt°-^g 

4.29 ± 0.18 ±0.20+0 2^ 



DGE 



Mx < 1.55 GeV 

Mx < 1.70 GeV 

P+ < 0.66 GeV 

Mx < 1.70 GeV, > 8 

Mx - g^ Pe > 1.0 GeV 

> 1.0 GeV 
p'e > 1.3 GeV 



35.3i|:3 



GeV" 



^4.8 
-4.8 

1+2.4 
-2.0 
7+3.5 
-3.2 

50.4l^3:g 



42.03 



24.43 
58.7j 



4.40 ± 0.16 ±0.12to.?9 
4.16 ± 0.18 ±0.14t^ J^ 
4.10 ± 0.19 ±0.17t^ J^ 
4.19 ± 0.19 ±0.12t;j;i^ 
4.40 ±0.16 ±0.181^13 
4.42 ±0.19 ±0.23+° " 
4.39 ±0.19 ±0.20 



0.13 
+0.15 



GGOU 



Mx < 1.55 GeV 

Mx < 1.70 GeV 

P+ < 0.66 GeV 

Mx < 1.70 GeV, > 8 GeV^ 

Mx ~ g^ Pf > 1.0 GeV 

p*i > 1.0 GeV 

p*e > 1.3 GeV 



-,+4.6 
-3.8 
j+4.2 
-3.6 
-1+8.6 
-6.3 
7+3.2 
-2.4 

60.213;° 
60.21^:° 

3+2.8 
^-2.3 



41.03 

46.81 
44.01 
24.71 



51. 



4.08 ± 0.15 ±0.1lt;|!;2? 
3.94 ±0.17 ±0.141:^;;;? 
3.75 ± 0.17 ±0.15l:o j^ 
4.17 ± 0.18 ±0.12t|^;25 
4.35 ± 0.16 ±0.181:;];?^ 



4.36 ±0.19 ±0.23 
4.33 ±0.18 ±0.20 



+0.09 
-0.10 
+0.10 



ADFR 



Mx < 1.55 GeV 

Mx < 1.70 GeV 

P+ < 0.66 GeV 

Mx < 1.70 GeV, g^ > 8 

Mx - g^ pI > 1.0 GeV 

p'l > 1.0 GeV 

pi > 1.3 GeV 



GeV" 



48.9i|^ 



62.01 
62.01 



-2.5 
-1+5.7 
-5.0 
-1+5.7 
-5.0 



3.81 ±0.14 ±0.11^0 20 

3.73 ± 0.16 ±0.13l:^;i^ 
3.56 ± 0.16 ±0.15t|];i^ 

3.74 ± 0.16 ±0.1llo:" 

4.29 ± 0.15 ±0.18+° }9 

4.30 ± 0.19 ±0.23to:w 
4.27 ± 0.18 ±0.19tg i^ 



Both limits arc larger than the theoretical limits, es- 
timated from D and Ds semileptonic decay rates, of 
3% [26, 27], and the more recent and stringent one of 
less than 2% [28, 29]. 

VII. CONCLUSIONS 

In summary, we have measured the branching fractions 
for inclusive charmless semileptonic B decays B — )- Xuii', 
in various overlapping regions of phase space, based on 
the full BaBAR data sample. The results are presented 
for the full sample, and also separately for charged and 
neutral B mesons. 

We have extracted the magnitude of the CKM clement 
\Vub\ based on several theoretical calculations. Measure- 
ments in different phase space regions are consistent for 
all sets of calculations, within their uncertainties. Corre- 
lations between \Vub\ measurements, including both ex- 
perimental and theoretical uncertainties are presented. 
They arc close to 100% for the theoretical input. 



We have obtained the most precise results from the 
analysis based on the two-dimensional fit to Mx~ q^, 
with no restriction other than p| > 1.0 GeV. The total 
uncertainty is about 7%, comparable in precision to the 
result recently presented by the Belle Collaboration [17] 
which uses a multivariate discriminant to reduce the com- 
binatorial background. The results presented here su- 
persede earlier BaBAR measurements based on a smaller 
tagged sample of events [11]. 

We have found no evidence for isospin violation; the 
difference between the partial branching fractions for 5° 
and B~ is consistent with zero. Based on this measure- 
ment, we place a limit on a potential contribution from 
Weak Annihilation of 19% of the total charmless semilep- 
tonic branching faction at 90% C.L., which is still larger 
than recent theoretical expectations [28, 29]. 

Improvements in these measurements will require 
larger tagged data samples recorded with improved de- 
tectors and much improved understanding of the simula- 
tion of semileptonic B decays, both background decays 
involving charm mesons as well as exclusive and inclusive 



TABLE IX: Summary of the fitted number of events A'„, the efficiencies, the partial branching fractions AB{B — >■ Xu£i^) and \Vub\ (10^^) based on four different QCD 
calculations of the hadronic matrix element as a function of the lower limit on the lepton momentum p^. The uncertainties on AB{B — ^ Xulv) are statistical and 
systematic, those for \Vub\ are statistical, systematic and theoretical. The uncertainties on all other parameters are statistical. \Vub\ values for BLNP and GGOU are 
not provided above 2.2 GeV due to large uncertainties. 
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decays contributing to the signal. Reductions in the theo- 
retical uncertainties are expected to come from improved 
QCD calculations for b — >■ uii/ and 6 — > S7 transitions, 
combined with improved information on the b quark mass 
and measurements of radiative B decays. 
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